#! /usr/bin/env python

import pyx
import pyfeyn
from pyfeyn.user import *

fd = FeynDiagram()
processOptions()

in1 = Point(-6,  2)
in2 = Point(-6, -2)
in11 = Point(-2,  2)
in22 = Point(-2, -2)

out = Point(4,  0)
out11 = Point(2, -3)
out22 = Point(2,  3)

vtx = Vertex(1, 0, mark=CIRCLE)

Fermion(in1,in11).addArrow().addLabel(r"\Pq")
Fermion(in2,in22).addArrow().addLabel(r"\Pq")
Photon(in11,vtx).addLabel(r"\PW")
Photon(in22,vtx).addLabel(r"\PW")
Fermion(in11,out22).addArrow().addLabel(r"\Pq")
Fermion(in22,out11).addArrow().addLabel(r"\Pq")
Higgs(vtx,out).addLabel(r"\PHiggs")


fd.draw("./VectorFusion_Higgs_production_Feynman.pdf")

